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Abstract. We present evolutionary synthesis models of starbursts on top of old stellar populations to investigate in 
detailed time evolution the relation between Ha luminosity and star formation rate (SFR). The models show that 
several effects have an impact on the ratio between L(Hq,) and SFR. Metallicity different from solar abundance, a 
time delay between star formation and maximum Hc-luminosity, and a varying stellar initial mass function give 
rise to strong variations in the ratio of Hq, luminosity to SFR and can cause large errors in the determination of 
the SFR when employing well-known calibrations. When studying star-bursting dwarf galaxies, and sub-galactic 
fragments at high redshift, which show SFR fluctuating on short timescales, these effects can add up to errors of 
two orders of magnitude compared with the calibrations. To accurately determine the true current SFR additional 
data in combination with models for the spectral energy distribution are needed. 
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1. Introduction 

The star formation rate (SFR) is one of the basic proper- 
ties of galaxies. It can be derived in different wavelength 
regimes (UV, optical, FIR) using empirical calibrations 
obtained from well studied samples of galaxy types, like 
spirals, irregulars, starbursts, or luminous IR galaxies. 
One of the most popular methods to derive SFRs from op- 
tical observations is by measuring the Balmcr line fluxes, 
which are very sensitive to the Hii regions surrounding 
massive young stars and therefore give a good measure of 
the very recent, shortlived, or ongoing SFR. When em- 
ploying Hq, one generally takes one of the simple linear 
calibrations as e.g. obtained by Hunter & Gallagher (1986, 
hereafter HG86) or [Kennicutt ct al. (1994| , KTC94) from 
observations and modeling of dwarf irregular and spiral 
galaxies, respectively. The calibrations are used in the 
forn^ 



L(Hq) [ergs-i] = • SFR [Mgyr-i] , 



(1) 



where the V^-factor is 1. 41 • 10'^'^ (HG86) or 1.26 • lO'^i 
(KTC94, [Kennicutt 199^ ) for a Salpeter IMF in the mass 
range of 0.1 to 100 Mq. 

These L(Hq)-SFR calibrations were derived for "nor- 
mal" galaxies with modest star formation rates. For sys- 
tems which behave differently from the ones for which 
the calibrations were derived, the application of this 
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^ The name of the variable = -£/(Hq ) /SFR is deduced from 
the fact, that the unit of actually is a squared velocity. 



method may not be appropriate. The validity of the SFR- 
determination in star-bursting galaxies, like e.g. blue com- 
pact dwarfs (BCDs), where these relations are frequently 
used, has never been shown. 

In fact, as the strongest output of ionizing photons 
is related to the most luminous, i.e. giant or supergiant 
phases of the ionizing stars, there may be a small time 
delay between abrupt changes in the SFR and the corre- 
sponding changes in the Hq flux. Leitherer et al. (1995[ ) 
showed that such a time dependence exists for stellar fea- 
tures in the UV. For systems with SF fluctuations on short 
timescales, the delay effect of the Balmer lines can also be 
important, e.g. for small scale systems like dwarf galaxies, 
where star burst durations are usually assumed to be of 
the order of a dynamical timescale, i.e. 10^ to 10^ years. 
SFR fluctuating strongly on short timescales may also 
have taken place in sub-galactic fragments before merg- 
ing together to hierarchically build up today's galaxies 
( plazebrook et al. 1999| ). 

We will show that indeed one may make very large er- 
rors when blindly applying the usual calibrations to small 
systems. We first present details on our model in Sect. |[ 
We then describe different effects that can affect the cal- 
ibrations for L(Hq) in terms of SFR, namely metallicity 
(Sect. H), short burst durations (Sect. and changes in 
the IMF (Sect. |). We finally summarize our results in 
Sect. |. 

2. Model description 

We use our evolutionary synthesis code specifically 
adapted to the modeling of starbursts in dwarf galaxies 
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(Kriiger et al. 1995; Weilbacher et al. 200C). It includes Table 1. Calibrations for different metallicities Z. 



specific modeling of gaseous emission lines and continuum 
based on the Lyman continuum photons emitted by hot 
young stars. 



We use the cu rrent Geneva stellar tracks (see Lejeune 
& Schi erer 2001 , for a recent compilation) in the metal- 
licity range from Z0/2O to 2^0. To be compati ble with 
HG86 and KTC94 we use the Salpeter IMF ( |Salpeter 
1955|yi)i the range of 0.15 to 85 or 120 Mq, as given by the 
tracks. We use the recently revised Lyman continuum pho- 
ton emission rates as given by Schaerer & de Koter (1997), 



who account for non-LTE effects, line blanketing, stellar 
winds, and the new temperature and gravity calibrations 
by Vacca et al. (199(:). To derive the luminosity, we 
sum up the emerging Lyman continuum photon emission 
7V(II") for each star, and convert it to a luminosity in Hq, 
using 



L(H„) [ergs^i] = 1.36 • lO"^^ 7V(H°) [s^^] 



(2) 



We have updated our models usin g the emission line r atios 
for low metallicities ob served by Izotov et al. (1997 ) and 
Izotov fc Thuan (1998 ) for a large sample of blue compact 
dwarfs (BCDs). 

We present two types of one-zone models. The first rep- 
resents a quiescent galaxy, where the SFR decreases slowly 
from the formation epoch with a timescale of 10 Gyrs. The 



other model was already discussed in detail by Weilbacher 
et al. ( pOOO ) in their interpretation of Tidal Dwarf Galaxy 



candidates, and c an also be used to model blue compact 
dwarfs (BCDs, see Kriiger et al. 1995 ). Here we put a star- 
burst on top of the stellar population of the undisturbed 
model. It is assumed to reach its maximum SFR of 20 
M0yr~-'^ after 10 Gyrs, and we vary the burst timescale 
tb from 10^ to 10^ yrs to model all the range of burst du- 
rations from small dwarf galaxies to mergers of giant gas- 
rich galaxies. We alternatively a bell shaped (Gaussian) 
burst or sharply rising and exponentially decreasing star- 
burst. We assume a minimum SFR after the starburst of 
O.lM0yr-i. 

3. Effects of Metallicity 

We successfully reproduce the calibrations of IIG86 and 
KTC94 using our quiescent models using solar metallicity 
and upper mass hmits of 120 and 85 M0, respectively. The 
values we derive for the four other metallicities are given 
in Tab. 0. 

While for solar metallicity models the agreement with 
IIG86's and KTC94's empirical calibrations is very good, 
it is also seen in Tab. |^ that both towards lower and higher 
metallicities the differences become significant. For subso- 
lar metal licities as e.g. in BCDs { {Z)^^^ = 0.002±0.001 = 



1/10^0, [Izotov fc Thuan 1998]) the difference in is 



found to be a factor of '^3.5. This is a result of low metal- 
licity stellar populations being both more luminous and 
hotter than a solar metallicity stellar population with the 
same IMF and mass limits. 
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Application of the empirical calibrations to estimate 
SFRs from H„ luminosities in low metallicity dwarf galax- 
ies hence can yield SFRs overestimated by a factor >3 due 
to metallicity effects. We also give in Tab. I the V^-value 
obtained in our solar metallicity model using the IMF of 
^calo (1986 ) instead of Salpeter's. Due to the smaller num- 
ber of high mass stars in case of a Scalo-IMF, the same 
SFR produces an luminosity that is lower by 35%. 

All other hydrogen lines, as e.g. Lq or Br.^, which can 
also be used to estimate SFRs, and their respective cal- 
ibration factors can easily be computed from our values 
for Ha and Eq. (|l|). 

At higher redshift, the [O ii]3727 line is frequently used 
to derive SFRs, as e.g. for Lyman Break galaxies. This 
line involves a direct metallicity dependence in addition 
to the one inherent in the hydrogen lines which is due to 
differences in temperature and luminosity of stellar pop- 
ulations at various metallicities. Since high redshift star- 
forming objects and, in particular, sub-galactic fragments 
both have low metallicities and possibly strongly fluctuat- 
ing SFRs, care is needed to derive SFRs from L([0 11]). We 
present the values of the Vp^j-factor in Col. 4 of Tab. |^ 
for the [Oii]3727 line, where the according linear relation 
with SFR is 



L([Oii]) [ergs-i] =Vfo„] - SFR [M 



(3) 



4. Short starbursts 



In Fig. |l| we show the time evolution of ^(Hq) for rapidly 
rising (instantaneous) starbursts with short timescales 
(tb = 10^ yr) for five metallicities. All models have their 
maximum SFR at a time of 10.0 Gyrs. We note two ef- 
fects: For lower metallicity Z the maximum H^ luminos- 
ity is higher than for high metallicity, e.g. by a factor of 
4.6 in case oi Z — 0.001 as compared with Z ~ 0.040. 
This is due to the higher temperatures of the low metal- 
licity stars, which more efficiently ionize the interstellar 
medium. It is also apparent that there is a delay of the 
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Fig. 1. L(Hq,) over time for instantaneous starbursts with Fig. 3. = L(Ha)/ SFR over time for Gaussian bursts 
short timescale of tb = 10^ yr. with different tiniescales and Z = 0.001. 
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Fig. 2. = ^(Hq)/ SFR over time for Gaussian bursts 
with short timescale of tb ~ 10^ yr. 

maximum in the luminosity with regard to the maxi- 
mum of the SFR. The offset between maximum SFR and 
maximum L(Hq) is higher for lower metallicities, 2.9Myrs 
for Z = 0.004 vs. l.SMyrs for Z = 0.040. This again is a 
result of the higher temperatures and hence the ionizing 
power of low metallicity stars. At low metallicity stars of 
lower mass and longer main sequence lifetimes contribute 
to the Lyman continuum emission. Since those take longer 
to reach their maximum Lyman continuum emission rates 
during their supergiant phase the delay becomes slightly 
longer than in the solar metallicity case. 

In Fig. H we plot the ratio between L(Hq,) and the 
SFR of our models as a function of time for Gauss-shaped 
bursts with a short timescale of tb = 10® yr. Before and 
after the starburst the empirical linear cahbrations agree 
well with the ratio seen in our solar metallicity model. 
With the onset of the burst the SFR rises faster than 
L(Hq) due to the delay in the maximum Lyman contin- 
uum photon production shown before in Fig. |l|. Therefore 
the V^-factor first decreases. After its minimum in- 
creases, because L(Hct) continues to rise until the ionizing 



stars have reached their supergiant phase. When the num- 
ber of Lyman continuum photons and therefore L(Hq) 
become maximal, the SFR has already decreased from its 
maximum value by a factor of '^3. After its maximum 
L(Hq) declines as the death rate of O stars is no longer 
compensated by SF. In this phase, however, the SFR de- 
creases even faster than the Hq, flux, and continues to 
increase. The sharp peak visible near 10.003 Gyrs, after 
which abruptly decreases, is a result of the constant 
minimum SFR of our models after the burst. If the SFR 
after the burst would go to zero, would diverge. The 
minimum SFR, which starts near 10.003 Gyrs, then acts to 
slowly bring down to values in agreement with those 
m Tab. H This time sequence during the starburst and 
the interplay between L(Hq) and SFR can also be seen in 
Fig. ^ below. 

It is obvious from Fig. ^ that for bursts with tb = 
10® yr strong discrepancies from the calibrations are seen 
for all metallicities, strongest at the lowest metal abun- 
dance, amounting to a difference of nearly two orders of 
magnitude as compared to the calibrations, even for solar 
metallicity models. 

In Fig. ^ we show the evolution of the V^-factor with 
time for five Gaussian starburst models with different 
burst timescales for our lowest metallicity Z = 0.001. 
With increasing burst duration, the maximum is shifted 
towards later times. It is strongest for bursts with tb = 
10® yr (which are shown in Fig. ^). This is a result of the 
convolution of the Gauss-shaped increase of the burst-SFR 
with the delay in maximum Hq emission due to the most 
massive and most luminous (= giant) stars. The decline 
after the sharply peaked maximum is an effect of the onset 
of the minimum SFR as discussed above for Fig. |[ The 
time evolution of significantly depends on the burst 
duration. For the longest burst with tb — 10^ yr the slow 
change in the SFR of the burst causes the delay effect to 
have negligible impact on V^. 

In conclusion we have shown that the SFR derived 
from Hq during or shortly after short bursts can be wrong 
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Fig. 4. L(Hq) over SFR for a bursting dwarf galaxy with 
metallicity Z = 0.001 and different IMF cuts. 



by factors ^10. Averaged over the entire burst duration, 
however, values for agree well with those from Tab. ^ 
for different metallicities. This applies to statistical anal- 
yses of samples of H ii galaxies. 

5. Effects of the IMF 

Kennicutt et al. (1994) have already investigated different 
V^-factors with different IMFs with identical mass cut- 
offs and found differences of an order of magnitude be- 
tween IMFs they used. Leitherer fc Heckman (1995 ) also 
presented properties like number of O stars, A^(H°), and 
EW(IIq,) etc. for instantaneous and continuous SF models 
with different IMFs. Here we want to show the effect of 
various IMFs on the time delay discussed above. 

Fig. ^ shows the luminosity £(IIq) plotted over the 
SFR of the model galaxy, while it experiences its starburst 
with a timescale of tb = 10^ yr. The calibrations of IIG86 
and KTC94 are plotted for reference. To indicate the evo- 
lution in time our "standard model" with Mup = 120 M0 
has additional dots each lO^yrs during the burst; the ar- 
rows show the direction of the loops in this diagram. 

It is obvious that the calibrations do not represent the 
real SFR very well. ^(Hq) at the peak SFR does fit very 
well with the calibrations for both models with Salpeter 



and Scalo IMF and a high mass cut at Mu 



120 M, 



0- 

But generally our models show that for an observed value 
of X(IIq,) there is an ambiguity between two values of 
the SFR. E.g. for L(Hq) = lO'^^ergs"^ the calibrations 
give a SFR of 7 ... 8 Mq yr~^, while one has to choose be- 
tween SFRs of and ~ 15 Mq yr'^ from our model with 
Salpeter IMF, and one needs additional data to determine 
the true current SFR. 

The evolution of the models after the starburst shows 
that for low SFRs (lower as one would expect from the 
calibrations) a galaxy could show high Hq, luminosity for 
quite some time (~ 3 • 10^ yr). For the two models with 
Salpeter IMF and lower mass cuts at Miow = 5Mq this 
effect is even more extreme. L(H„) at the peak SFR is 



already underestimated by the calibrations by factors of 3 
and 5, respectively, for high mass cutoffs of Mup — 40 and 
120 Mq. For a given Ha-luminosity the calibrations yield 
a SFR too high by at least one order of magnitude when 
compared to the values from our models. 

6. Conclusions 

When observing small scale star-forming entities like 
dwarf galaxies or sub-galactic fragments, where dynam- 
ical timescales and hence burst durations may typically 
be of the order of 10® yrs or less, one should be aware that 
the currently used calibrations on the basis of Hq lumi- 
nosities may yield SFRs with large errors. These calibra- 
tions were determined from larger low-level star-forming 
systems (where they work very well), and when applying 
them to star-bursting dwarf galaxies, errors of factors 3 to 
100 may affect this determination of the SFR. The ratio 
of Ha-luminosity to SFR depends on the metallicity of the 
object, and the age and the duration of the starburst. 

To more accurately determine the true current SFR 
of a (star-bursting) dwarf galaxy additional information 
about the spectral energy distribution (SED) is needed. 
Observations in at least three optical/NIR filters or a spec- 
trum with sufficient wavelength coverage to determine the 
slope of the SED could be used in comparison with models 
that include gaseous emission to eliminate ambiguities in 
the relation of the Hq luminosity to the SFR. 
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